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Chapitre 1 

Introduction 



C'est la troisieme fois que vous allez etudier Pelectrodynamique. En effet, vous 
avez deja pu vous familiariser avec ces concepts au gymnase ainsi que lors du 
troisieme semestre de physique generale. Quel est done l'objectif de ce cours ? 
Bien sur, vous retrouverez ici quelques rappels des notions deja acquises. Mais ce 
sera egalement l'occasion d'aborder des sujets qui n'ont pas encore ete traites. 
D 'autre part, ce cours n'a pas la pretention de couvrir tous les aspects de Pelectrodynamique. 
Nous nous limiterons aux deux problcmes suivants : 

• Pelectrodynamique classique comme base de la relativite restreinte, 

• approche microscopique des phenomenes electrodynamiques macroscopiques. 

Nous vous invitons egalement a consulter les references suivantes : 

• Classical electrodynamics, J.D. Jackson, John Wiley & Sons 

• The Classical theory of fields, L.D. Landau & E.M. Lifshitz, Pergamon Press 

• Electromagnetisme, R.P. Feynman, InterEditions 
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Chapitre 2 

Solution des equations de 
Maxwell en presence de 
charges et de courants : 
fonction de Green 



2.1 Les equations de Maxwell, le systeme d'unites 

Tout phenomene d'electrodynamique classique est decrit par les equations de 
Maxwell : 

V-£=- (2.1) 

VxE = ~ (2.2) 
V-B = (2.3) 

C 2 V x B = - + (2.4) 
eo ot 

ou E(x) est le champ electrique, 

B(x) est le champ magnetique, 

p{x) est la densite de charge, 

j(x) est la densite de courant, 

c est la vitesse de la lumiere, c = 299 792 458 m/s, 

e est la permittivite du vide, e = = 8.85 • 10" 12 ^ 

A ces equations, nous devons egalement ajouter la force de Lorentz : 

F = e{E + vxB) (2.5) 

oil e est la charge electrique de la particule. 

Toutes ces equations sont ecrites dans le systeme international d'unites MKS 
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(Metre Kilo Seconde). L'unite de courant en est Tampere A, et l'unite de temperature 
est le Kelvin K. 

• Definition de l'ampere : un courant est de 1 A si la force agissant entre 
deux lignes electriques paralleles parcourues par le meme courant est de 
2 • 1(T 7 N/m. 

• L'unite de charge est le Coulomb C : 1 C = 1 A ■ s. 

• L'unite du potentiel electrique est le Volt V : V = ^ 

• L'unite du champ electrique : V/m. 

• L'unite du flux magnetique est le Weber Wb : Wb = V ■ s. 

• L'unite de la densite de flux magnetique est le Tesla T : T = ^5 = ^# 
(1 T = 10 4 gauss). 

• L'unite de la resistance est le Ohm 17 : 17 = ^. 

• L'unite de capacite est le Farad F : F = y 

2.2 Les limites de l'electrodynamique classique 

Comme toute theorie physique, l'electrodynamique classique a ses limites d'ap- 
plications : elle n'est pas valable a tres petites distances et a tres hauts champs, 
a cause de la physique quantique. Elle pourrait egalement ne plus etre valable a 
tres grandes distances. 

De l'electrodynamique, nous savons que toute charge en mouvement rayonne. Ceci 
nous conduit a la conclusion suivante : un modele planetaire pour l'atome est 
instable. II faut done batir une theorie quantique de l'atome pour resoudre ce 
probleme. A ce sujet, nous vous renvoyons aux series d'exercices. 
D'autre part, nous ne savons pas experimentalement si l'electrodynamique clas- 
sique reste correcte a de tres grandes distances, / > 10 19 m. 

II n'en reste pas moins que l'electrodynamique classique a une grande region d'ap- 
plicabilite. 

2.3 Rappels des cours passes 

Dans ce qui suit, nous allons considerer que vous connaissez toutes les operations 
d 'analyse vectorielle. Nous rappelons : 
• en cordonnees cartesiennes : 
gradient : 




divergence : 
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rotationnel : 

V x A= - ^] ex + - ^] e 2 +(— - — ] < ■ 

\dx2 dx% J \dxs dxi J \ dx\ dx^ 

D 'autre part, nous considerons que vous connaissez : 
• Le theoreme de Gauss : 

J V • A dV = i . \ ■ liS 



Le theoreme de Stokes : 



j> A-dl = j (V x I) ■ dS 



De cette facon, nous pouvons recrire les equations de Maxwell sous forme integrale : 

E ■ dS = 9- 

"Le flux de E a travers une surface fermee est egal a la charge a l'interieur de 
cette surface divisee par eo" . 

E-A = ~ I B-dS 
r dt 

qui est la loi de Faraday, 

' B ■ d & = 



qui traduit la non-existence de monopoles magnetiques, 

c 2 j> B-dl = ]r f ^'^ S+ Wt f E- dS 

V v ' 

loi d' Ampere 

Dans les sections qui suivent, nous rappelons les diverses methodes que vous avez 
deja vues pour resoudre des problemes simples. 

2.3.1 Methode pour l'electrostatique 

Nous avons le potentiel eft qui est defini par E = qui doit resoudre 

V x E = 0. Le potentiel verifie l'equation de Poisson : 

A0 = — ^ (2.6) 



est defini a une constante pres, une solution de l'equation [2.6] est : 

1 

4vre 



m = ^- 1 p -^i (2.7) 
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2.3.2 Methode pour la magnetostatique 



Nous avons le potentiel vecteur A qui est defini par B = V x A, qui doit 
resoudre V • B = 0. Nous voyons done que le potentiel vecteur A n'est donne qu'a 
un gradient pres, puisque A + Va donne le meme champ magnetique. II nous faut 
done nous donner une jauge, prenons par exemple la jauge de Coulomb : V- A = 0. 
Les equations de la magnetostatique ont egalement la forme de l'equation de Pois- 
son 

AA = — ^ (2.8) 
ces equations ont une solution du type 

4ire c A J \x-x'\ 

Les equations [2.7] et [2.9] donnent une solution au probleme suivant : nous avons 
une distribution connue de courants et de charges, trouver A et <f>. 



2.3.3 Phenomenes dependant du temps : conservation de la charge 

Une consequence importante des equations de Maxwell est la conservation de 
la charge electrique. En effet, en utilisant l'equation [2.1], nous avons 



D' autre part en prenant [2.4] : 



c 2 V'(VxB) = ^|-V'E (2.11) 
v / e ot 



=o 

En combinant [2.10] et [2.11] nous obtenons : 

dp 



. )t +V-j = (2.12) 



qui est appelee equation de continuite. 



2.3.4 Methode pour l'electrodynamique 

II s'agit maintenant de resoudre les equations [2.2] et [2.3] simultanement. Nous 
avons vu que nous avons une certaine liberte dans le choix des potentiels (a une 
constante pres). II faut done que les solutions soient invariantes de jauge : 

A — > A + Va 
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v v dt 
Nous allons verifier cette invariance : 

B' = Vx(I+ Va) = Vxl + VxVa = B = Vxl 

=o 

-_(^ + Vo)-V^-^j= £ =-^-V* 
<=<"' (V| - »V) a = 0. 

Une remarque vient alors, nous pouvons choisir la jauge de fagon a simplifier les 
equations de Maxwell. La jauge de Lorenz 1 donne, par exemple, des equations 
decouplers. Le choix est le suivant : 

Nous obtenons les equations suivantes pour les potentiels : 

"^+4S = ^ (2.14) 
c 2 ot 2 eo 

-Al+^^ = ^2 ( 2 - 15 ) 
c 2 ot 2 eoc 2 

Si la densite de courants et la densite de charges sont nulles, les equations [2.14] 
et [2.15] sont les equations d'ondes avec la vitesse c. 

Le prochain objectif est de trouver des solutions pour ces equations en utilisant la 
methode des fonctions de Green. 



2.4 Methode des fonctions de Green pour les equations 
different ielles 

Dans ce paragraphe, nous allons considerer les deux equations suivantes : 

_A0=^ (2.16) 

qui correspond au cas statique, et 

ia 2 ^ P (x,t) 

-A0+^--2- = (2.17) 

c 2 ot 2 eo 



1 Rem : Cette jauge est nommee d'apres le physicien danois Ludvig Lorenz (1829-1891), alors 
que le physicien neerlandais Hendrik Antoon Lorentz (1853-1928) a donne son som aux transfor- 
mations de Lorentz et a la force de Lorentz. 
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pour le cas dynamique. 

Definition : Une fonction de Green est une fonction de (x, x') [ou bien de (x, t; x', t')\ 
qui satisfait les equations : 

- A g G(x, a?) = 5 3 {x - a?) (2.18) 

dans le cas statique, et 

(~ A *+^^) G(x,t;x\t') = S 3 (x-x')6(t-t') (2.19) 
dans le cas dynamique. 

Vous remarquez l'apparition, dans ces deux dernieres equations, d'un nouvel objet : 
la fonction 5 de Dirac. Ce n'est pas a proprement parler une fonction, mais une 
distribution. 

L'idee qui se cache derriere ce nouveau concept est la suivante : comment pouvons- 

nous decrire une charge ponctuelle q ? 

Une possibility est d'utiliser la fonction 5 de Dirac : 

avec v = fvre 3 , ainsi nous avons p e (x — x') = q5 3 (x — x'), et plus formellement : 

lim p e (x — x') = q5 3 (x — x') 

La forme explicite de p e n'est pas essentielle, plusieurs representations de la fonc- 
tion 5 peuvent etre prises, par exemple : 

S 3 e (x - x>) = (±) 2 exp (Jl^l) (2 . 21) 



7re z \ e 



La fonction 5 peut etre definie dans des espaces de toutes dimensions, nous ecrirons 
pour un espace a n dimensions : S n (x — x'). 



2.4.1 Proprietes de la fonction delta 

Nous listons ici quelques proprietes de la fonction 5 : 

a) 

h = J f(x')6 3 (x-x')d 3 x' = /(£) (2.22) 

• Preuve : 
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h = lim / f(x')5't(x - x')d 6 x' = lim / f{x - z)5 3 {z)d 3 z (2.23) 

e^O ./ e^O . 



ou nous avons fait le changement de variable x — x' = z. Ce qui donne en 
developpant : 

h = lim \ (fix) -z-Vf + ..) -d 3 z 



( 



I\ = lim 



v J\g\ <e v 



zd 3 z+... 




b) 



h = J f(x')V-,5 3 (x-x')d 3 x' = -V/(x) 



Preuve : 



h = lim / fix')V -Mix- x')d 6 x' 
e->oJ x 

En utilisant la formule du gradient et a) : 

I 2 = I f(x')6 3 {x - x')d~S -V/(x) 



car 5 3 (x — x') = a une surface infinie. 
D'autres proprietes seront demontrees en exercice. 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



Revenons maintenant aux fonctions de Green. Pourquoi sont-elles si utiles ? Si 
G est connue, alors, nous avons tout de suite une solution pour les equations 
[2.16] et [2.17] : 

En effet, pour l'equation [2.16], 



(x) = — I G{x,x') P (x')d A x' 
eo 



appliquons l'operateur —A sur l'equation [2.29] : 



A0(s) 



1 

eo 



-AgG(x, x')]p{x')d 3 x' 



(2.29) 



(2.30) 
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qui donne par definition de la fonction de Green 

1 



A(f)(x) 



eo 



8 z {x-x')p{x')d 6 x' 



(2.31) 



L 'equation de Poisson est done bien verifiee par la solution ainsi trouvee. 

Nous pouvons de la meme facon ecrire une solution pour le cas dynamique (pour 

l'equation [2.17]) : 



(x,t) = — / G(x,t;x',t')p(x',t')d 3 x'dt' 
eo 



(2.32) 



Comme plus haut, nous pouvons appliquer l'operateur de d'Alembert pour 
retrouver, cette fois-ci, l'equation [2.17] : 



i d 2 ' 




c 2 dt 2 


eo 


i d 2 '- 




c 2 dt 2 _ 


eo 



2 1 



1A_ 



G(x, t; x', t')d d x'dt' (2.33) 



5 3 (x - x')<5(* - t')p(x', t')d 3 x'dt' = -p(x, t) (2.34) 



eo 



Comme vous remarquez, il est plus facile de trouver une fonction de Green que 
de resoudre une equation generale. En effet, dans notre cas, la distribution a une 
forme tres specifique : elle est egale a zero sauf dans un point. II faut egalement 
noter que la fonction de Green n'est pas definie de facon unique : nous pouvons 
toujours lui ajouter une solution de l'equation homogene. 

2.4.2 Exemple 1 : cas statique 

Pour faire un premier exemple, nous allons chercher la fonction de Green pour le 
potentiel scalaire, dans le cas de l'electrostatique. Nous devons done verifier 
l'equation suivante : 

- A £ G(x, x') = 5 3 (x - x 1 ) (2.35) 

pour une charge ponctuelle situee en x' = 0. Nous considerons la solution de 
symetrie spherique ; nous rappelons la forme du Laplacien en coordonnees 
spheriques 

i d 2 



A 



1 d ( 9 d 



+ 



1 d 



d 



sin ( 



+ 



r 2 dr \ dr J r 2 sin 9 88 \ 36 J r 2 sin 2 6 dip 2 
Si nous prenons la representation S e , alors l'equation [2.35] devient pour r > e : 



1 d ( ,dG 



o 



(2.36) 



en simplifiant, nous obtenons 

2 dG 



_ dG Const 
r"— = Const => — = r— 
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Finalement 



G(r) 



C 



\r>e 



(ou C est une constante differente de Const.) 

II nous reste a determiner cette constante C, grace a l'equation [2.35] : 



.1JL ( r 2dG) A _* 
r 2 dr \ dr J 



4ir 



4tt 
/ 

V 



d_( 2 dG 

dr I dr 



! l.nr'dr — 1 

2 dG 



C 




-r 



dr 



AttC = 1 



Nous avons done C = -h-. Nous en deduisons la forme de la fonction de Green : 



G(x, x') = 

4ir\x — x'\ 

Ainsi, nous pouvons en tirer la solution pour l'electrostatique : 



(2.37) 



1 

4vre 



p(x') 



1 



\x — x 



-d 6 x' 



(2.38) 



2.4.3 Exemple 2 : cas dynamique 

Nous pouvons effectuer le meme genre de calculs pour le cas dynamique. 
L'equation a resoudre devient : 

G{x, t; x', t') = 5 3 {x - x')8(t - t') (2.39) 

Nous allons, pour resoudre ce probleme, faire des transformations de Fourier. 
Nous les definissons de la facon suivante : 

i r+oo 

G(x,t;x',t') = — G(x,uj;x',t')e- iujt duj (2.40) 
ce qui donne pour la transformation inverse : 

/+oo 
G(x,t;x',t')e iujt dt (2.41) 
-oo 



-A, 



1A 
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Pourquoi passons-nous dans l'espace de Fourier? L'interet en est evident lorsque 
nous ecrivons les relations de transformations [2.40] et [2.41] pour la fonction 5 : 



S(t - 1') 



1 

2^ 



-Mt-t') duJ 



car 



e luJt 8(t)dt = 1 



En reportant ce resultat dans l'equation [2.39] nous arrivons a : 



-A, 



G(x,uj;x',t') = 5 3 (x-x')e iujt ' 



(2.42) 
(2.43) 

(2.44) 



Posons K = — et prenons x' = en utilisant la symetrie spherique. L'equation 
[2.44] devient alors : 



Si r > e alors [2.45] devient : 



5 3 e iutf 







(2.45) 



(2.46) 



pour resoudre [2.46], ecrivons la fonction de Green sous la forme : G = ^F(r). 
Nous avons : 



dG 
dr 



^ + — -» I" ( r2 ir) = - F ' + F ' + rF " = rF " ( 2 - 47 ) 

r z r or \ or J 



En injectant cela dans [2.46], nous obtenons l'equation differentielle bien 
connue : \ (F" + K 2 F) = 0, dont les solutions s'ecrivent F = Ae iKr + Be~ iKr . 
Nous avons done la forme generale de la fonction de Green : 



G = -e iKr + -e~ iKr 



(2.48) 



• Considerons le cas r < e : 
comme nous avons maintenant : 



[ K 2 Gd 3 x ~ 0(K 2 )-e 3 ~ 0((Ke 
Jo e 



) 2 )^0 (6-0) 



nous trouvons 



A + B = e 



, 1 

47T 



la solution B = est appelee fonction de Green retardee, notee Gr, 
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la solution A = est appelee fonction de Green avancee, notee Ga- 

Nous verrons plus tard le pourquoi de ces appellations. Prenons pour l'instant la 

solution retardee : 

1 /" +0 ° e iK\x-x'\+iu>(t'-t) J 

G R (x,t;,x',t') = ±- / - (2.49) 

^ J-oo \X-X'\ 47T 

G R (x,t;,x',t') = ! +0 ° e^-^(t'-t)^ (2 . 50) 

2vr \x - x '\ J-oo 47r 

Par definition de la transformee de Fourier de la fonction 5, nous trouvons : 

G R (x,t;,x',if) = — 1 —8 ( t'-t + 1 \x-x'\ ) — (2.51) 
\x - x'\ V c / 47r 

Par pure analogie, nous pouvons ecrire la solution avancee : 

G A {x,t;,x',t') = l -^5\t' - t--\x- x'\ ) — (2.52) 

\x - x'\ \ c J 47r 



Ainsi, la solution pour le potentiel electrodynamique s'ecrit : 

A(x ,t) = - / p(x',t') — -^5 ( t' - t + -\x - x'\ ) (fx'dt' (2.53) 

47re y v \x-x'\ \ c ) 

cj)(x,t) = - / p{x',t- -\x-x'\) —d 3 x' (2.54) 

4vre J c |f _ x '\ 

Ce qui signifie que le potentiel au temps t depend de la distribution de charges a 

l'instant t — \\x — x'\. Nous avons done un signal retarde. 

Pour le potentiel vecteur, la solution peut-etre trouvee immediatement : 




(2.55) 



Remarque : ces solutions sont des solutions particulieres. Une solution generale 
est une solution particuliere + une solution arbitraire de l'equation homogene. 
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2.5 La radiation du dipole (rappel) 



Une consequence importante des solutions retardees est l'emission d'ondes 
electromagnetiques. Le resultat principal est la radiation du dipole. 




Nous considerons le cas suivant : r <C R et v <C c, ou v represente la vitesse des 
particules dans Pemetteur. Les potentiels sont alors donnes par les formules : 



A(R,t) 



P\t-« 



47re c 2 R 



(f>{R,t) 



1 



R 



4tt€ R 3 
oil p est le moment dipolaire du systeme, qui s'ecrit : 



Prenons l'exemple de deux particules ponctuelles de charges opposees 



-e 



+e 



(2.56) 



(2.57) 



(2.58) 



alors la densite de charges est p = —e5 3 (x + |) + e<5 3 (x — j), nous avons done 
pour le moment dipolaire : 

d d r 

P = e- + e- = ed 

Pour des distances tres grandes, les champs s'ecrivent alors : 

1 lfixR 1 



B = 



4:TT€qC 2 c R? 



0( 



R J 



-* 1 (px R)x R - R 

E = -, o— — ^5 = cB x — 

4ire c 2 R 3 R 



(2.59) 
(2.60) 
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Le vecteur de Poynting se calcule comme suit : 

S = e c 2 E x B 
le flux d'energie est done donne par : 

— — = f S • da = 5- ( 

at J otreocr 
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Chapitre 3 



Les potentiels de 
Lienard-Wiechert 



Considerons le probleme suivant : 

Une particule chargee suit une trajectoire arbitraire, nous voulons connaitre les 
champs crees par cette particule, l'observateur etant place a l'origine. 



Notons ici que ce probleme est a l'origine de : 

- la construction de la mecanique relativiste, 

- la comprehension de tout probleme de radiation, pour un cas arbitraire (v ~ c). 

3.1 Calcul des potentiels 

Nous prendrons les notations suivantes pour la trajectoire de la particule : 



z 




X 




(3.1) 
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Soit, pour le vecteur position R(t) = R2, R3). Avec ces notations, nous 
pouvons ecrire la densite de charges et la densite de courants comme : 



p = eS(x-R(t)) 



(3.2) 



j = eR5(x - R(t)) = ev5(x - R(t)) 



(3.3) 



II est maintenant aise d'ecrire la forme des potentiels. Nous avons, par exemple, 
pour le potentiel scalaire : 



<f>(0,t) 



— — — I eb \ a? 
4vreo J \ 



R(t 




d 3 x' 



(3.4) 



expression compliqucc 



L'idee qui nous vient est d'effectuer un changement de variable, de fagon a 
simplifier l'expression du potentiel. Faisons ce changement en posant : 



x = x' -R(t- ±) 



(3.5) 



Nous avons done : 



d 3 f=^5-d 3 s' = det 
D(x') 



dx 



dx' 



d 3 x' 



(3.6) 



ou nous n'oublions pas le Jacobien de la transformation. Nous utiliserons le fait 
que ^r\x\ = p. Faisons le calcul explicitement : 



D(x) 

— = det 
D(x') 




dxi 
dx' 2 
8x2 
dx' 2 
8x3 
dxi, 



dxi \ 
dx' z \ 
8x2 
dx' 3 
8x3 
dx'.. 



(3.7) 



D(x') 



Vl 



c\x'\ 



V2 
V3 



c\x'\ 

A 

c\x'\ 



Vi 



c\x'\ 



1 + W2 T% 

c\x I 



^3 



c\x'\ 



Vl 



V2 



c\x'\ 

x 3 

c\x'\ 



det 



ViX i 
hi + -t=t 
car 



(3.8) 



Cette derniere equation peut etre evaluee plus facilement dans un systeme de 
coordonnees pour lequel nous avons : Ri = R2 = 0, -R3 7^ 0. Nous pouvons, en 
effet, toujours nous ramener a ce point par une rotation ; nous obtenons alors : 



1 
1 








l + v *W\ 



1+V3-4,- 



1 _ X' 

1 + v ■ — 



c\x' 



c\x 



(3.9) 



nouveau systeme systeme arbitrairc 
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Ainsi, nous pouvons ecrire explicitement le potentiel scalaire, apres le 
changement de variable [3.5]. Nous injectons l'expression du Jacobien [3.9] dans 
[3.4], et nous obtenons : 



0(0, t) 



1 



47rer 



S(x) 



Ml i + 

1 ' c\x'\ 



-d 3 x 



(3.10) 



ou x' est solution de l'equation [3.5]. Done : 

j . . 1 e 1 e 
0(0,*) 



1 



47re i + f^. Ml 47re M| + -u ■ x' 



(3.11) 



avec x' = fl(t_ JfLl) = i?(t'). Nous avons alors la forme finale du potentiel 
scalaire : 

0(O,t) = -^ \ ^_ (3.12) 

47re 0j R(t') + ^ ■ R(f) 



Pour le potentiel vecteur, nous obtenons : 

- e? 1 



47re c 2 i?+ Itf. ij 



(3.13) 



t'=t- 



Ces expressions traduisent, du point de vue du sens physique, le retard. Les 
phenomenes electrodynamiques ont un temps de propagation fini : . Pour 
nous convaincre de la validite des formules trouvees, choisissons une vitesse de 
propagation infinie (c — ► oo). Nous retrouvons alors la loi de Coulomb pour 0, et 
l'expression du potentiel magnetostatique pour A (loi de Biot-Savart). 

3.1.1 Derivation des champs magnetiques et electriques pour le 
cas non-relativiste 

Nous allons partir des expressions des potentiels scalaires et vecteurs : 

4>(x,t) 



4?re o R+^v-R 



(3.14) 



ou nous avons : R = Ro — x, avec Ro(t) definissant la trajectoire de la particule 
chargee. Nous trouvons alors t' grace a l'equation suivante : 



t' = t 



R(t>) 



(3.15) 



Nous voulons trouver E et B dans le cas non relativiste (- <C 1), nous allons 
done negliger les termes du type : (^) 2 , (^) 3 , vv, etc. D'autre part, nous ne 
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considererons que les champs loin de la particule, nous negligerons done 
egalement les termes du type (-^) avec n > 2. 

- l ere etape : solution de l'equation [3.15] 

une premiere approximation est de dire t' = t, 

dans un second temps nous ecrivons t' = t — = to(t, x), 

enfin nous pouvons ecrire : t' = t — — - c - , qui se developpe comme suit : 

i! Pztt— + ^i?^— = to + 0(v), nous n'aurons pas besoin des termes en 
0(v) pour t', voir plus bas. 

- 2 eme etape : approximations pour R et v 

R{t') = R(t ) + R ■ 0{v) = R(t ) + 0(v 2 ) 



etape : developpement de <f> et A en puissances de v. 

to (3.16) 



^«5^S(1 + 0(«))U 



4 eme etape : calculs intermediaries 

comme E = — V(j) - |^ et B = V x i, nous avons besoin de calculer et 



9x 3 % 8 s, -V c^-tj-M (3JT) 

Ui + 0(v a ) (3.18) 



<9t> i dvi dto . , / w L . , 



etape : Notons que : 



de meme : 



-(-)=o(±- 

dx\R \R 2 



-(-)=o(± 

dt\R \R 2 



done ces termes peuvent etre negliges a grande distance. 
Derniere etape : calcul de E et B 
pour le champ magnetique : 

B =e d ^=e- 6 1 dVh 
1 u dxj ^ 4ireoc 2 R dxj 
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ce qui donne, en utilisant le calcul [3.17] : 



B = -—^-^ifxR (3.19) 
47re c 3 R 2 



pour le champ electrique : 



E _ d4> dAj ^ e 1 1 ( dv g\ e i)i 



dxi dt 4ireo c R 2 \dxi J 4ireo c 2 R 
de la meme fagon, nous utilisons [3.17 et 3.18], et nous obtenons ainsi : 

dans ces expressions le temps est to. Rappelons que ces expressions ne sont 
valides que pour R — > oo. 
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3.1.2 Comparaison avec le dipole 

Dans cette section, nous allons chercher les explicites expressions 
non-relativistes, pour les potentiels <fi et A. Nous les comparerons alors aux 
expressions obtenues pour le rayonnement du dipole. 
Cherchons une solution par iteration de l'equation 



t' = t 



R(t>) 



(3.21) 



de laquelle une premiere approximation est : t' = t, 

en seconde approximation, portons l'expression t' = t — ^r^, 



dans 



R-v 



R{t') + 



R, R-v 

R{t--) + 

c c 



(3.22) 



(3.23) 



en faisant un developpement de Taylor au premier ordre en ^, nous avons pour 
[3.23] 



R(t) 



dRR R-v 

dt c c 



Faisons ici un petit calcul intermediate : 



dt dt 



2R 



l n 8R = R-v 
2R dt R ' 
Grace a ce resultat, nous pouvons ecrire [3.24] : 

. R ■ v R R-v „ . . 

jK«)-_ 7 + _ = jK t) 

ce qui donne pour le potentiel scalaire : 

e 1 



4vre R{t) 



(3.24) 

(3.25) 
(3.26) 

(3.27) 
(3.28) 



et pour le potentiel vecteur : 



A 



v \t-* 



4vre c 2 R(t) 

II nous reste a retrouver les expressions du dipole rappelees en [2.56] et [2.57]. 



(3.29) 
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Pour cela, nous considerons deux charges +e et — e, situees en R\ et i?2- Les 
deux charges sont proches l'une de l'autre, de telle sorte que si R = \{R\ + R2), 
nous ayons |.Ri — -R2I «A En utilisant le principe de superposition, nous 
trouvons d'apres les equations [3.28] et [3.29] : 

1 P\t-S 

- 1 (3.31) 



r 4ire R 3 (t) 
avec la definition du moment dipolaire : p = eR\ — ei?2- 

L'equation pour A correspond exactement a l'equation dipolaire. Prouvons qu'il 
en est de meme pour <p. 
Nous remarquons que : 



p + — p 



p(t)-^~ + ~P=P(t) (3-32) 



L'equation [3.31] est done la meme que [2.57]. 

Ces transformations sont importantes, car elles nous permettent de deduire les 
champs des equations de Maxwell : 

dA 

E = -V0 - (3.33) 
B = V x A (3.34) 
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3.2 Cas simple : particule le long d'une ligne droite 

Pour mieux comprendre le sens physique des potentiels de Lienard-Wiechert, 
nous allons etudier le cas d'une particule en mouvement le long d'une ligne 
droite. 

Regardons le schema suivant : 




La trajectoire est parametrisee par : 

' x(t) = 

< y(t) = o 

( z(t) + 

et nous nous trouvons en un point arbitraire x. 
Ecrivons les potentiels dans ce cas 



4> = 



l 



4^eo y/x 2 + y 2 + {z- z{t')) 2 + ±z(z(t') - z) 

" v ' 

=R 



(3.35) 



A, 



z 
7i 



tout etant evalue a l'instant 



1 



t' = t- -^x 2 + y 2 + {z-z{t')) 2 



(3.36) 



(3.37) 



Cette derniere equation est une expression compliquee de t' en fonction de 
dependances arbitraires en R et t. Essayons maintenant de calculer les champs 
magnetiques et electriques. 



3.2.1 Champ magnetique 

Le champ B est donne par B = V x A. Dans notre hypothese, seule la 
composante A z est non-nulle, done 




B= \ - d -t \ (3-38) 
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Nous avons une symetrie cylindrique, done seule la composante du champ 
magnetique est non-nulle. II est done suffisant de prendre y = et de calculer 

D v — dx ■ 



3.2.2 Champ electrique 

Le champ E est donne par E = —V<p 
composante A z est non-nulle, done 



Dans notre hypothese, seule la 



E 



( 



(3.39) 



d± _ dA z 
dz dt 



Nous avons une symetrie cylindrique, comme plus haut, il est suffisant de 



prendre y = et de calculer E x 



dx 



et E 7 



_ _ dAz 
dz dt 



Nous arrivons ici au coeur du calcul : t' est une fonction de x, z, t, si nous 

prenons y = dans l'equation [3.37]. 

df_ df_ 8f_ 

dx ' dz ' dt 



Nous devons calculer :— — — 



3.2.3 Calculs intermediaires 

Nous allons dans cette partie evaluer les differentielles. 



dx 



Done : 



x z — z(t') , ..dt' 
R + R [ ~ Z) dx~ 



d£_ 
dx 



lx_ 
c R 



1 + 



R-v 

cR 



X R - V dt' 

R + ~R~dx 



(3.40) 



(3.41) 



De la meme fagon : 



Ce qui donne : 



d£_ 
dz 



z - z(t') z- z(t') 



R 



R 



1^ 

dz 



dt , _i_z_-4t> 1 



R 



dz 



1 + 



R-v 

cR 



(3.42) 



(3.43) 



II nous reste a calculer : 

dt' _ 
~dt ~ 



lv-Rdt' 
c R ~dt 



(3.44) 
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Soit : 



dt' 
~dt 



1 + 



1 v-R 
c R 



(3.45) 



3.2.4 Resultats 

Les equations [3.41, 3.43, 3.45] peuvent etre injectees dans les expressions des 
champs [3.38, 3.39], nous trouverons ainsi les resultats cherches. 
Calculons la composante E x 



E ,■ 



1 



4vre I R+ vji 



R 



+ 



R c c ax 



4vre f R + ti-R x 2 



' cR 



x I v ■ R v 2 v ■ R\ 



tvre ( R + ^ J R^ 2 R 



1 1 



1 + 



Nous obtenons finalement : 



E ,■ 



( v ■ R v 2 v-R 
vjl \ R c c 

cR \ / 



v 2 v-R 



4:71-6, 



°(R+^y 



c 2 c 2 



Ey=0 



(3.46) 



(3.47) 



Pour la composante E z , il suffit de faire un calcul similaire. Le resultat final, 
pour le cas general (trajectoire arbitraire) est le suivant : 



E 



47T60 



1 



R+-R 



R x 



R+ZR) x v 



(r+*£Y v c j c 2 (r + *- r 

Pour trouver le champ magnetique il suffit de faire : 

1 R 



B 



c R 



x E . 



(3.48) 



(3.49) 



3.2.5 Interpretation des resultats : la radiation 

Nous allons maintenant regarder le sens physique des equations [3.48, 3.49] que 
nous venons d'obtenir. Pour cela, nous allons etudier la limite lorsque R — > oo. 
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Dans cette limite, nous distinguons deux termes dans l'expression [3.48], ils ont 
deux comportement differents 



E 



O 



l er terme 



1 
i?2 



exactement comme pour une charge ponctuelle. Alors que le second terme 
donne : 

2 e terme 

terme qui est dominant aux grandes distances. Le champ B qui est lie au champ 
electrique a le meme comportement, ce qui signifie que le vecteur de Poynting se 
comporte en O (^)- 

Done le flux d'energie « R 2 ■ O (jp) ~ const. Alors, il y a de la radiation 
d'energie. 

Remarque : le second terme est different de zero seulement si v ^ 0, e'est-a-dire 
si la particule est acceleree. 



Calculons le vecteur de Poynting pour des distances R importantes : nous avons 
pour le champ electrique, dans le cas general : 



E 



R x 



R+^R) x v 



R-^oo 47re ( 



o 



cHR+^f 



(3.50) 



Ce qui nous dit que E _L R done E ■ R = 0. Or, le vecteur de Poynting est donne 
par 



S = e c z E x B : 
Soit, en recrivant le double produit vectoriel : 



1 1 



e c---E x [RxE 
R c 



S 



R 



R(E) — E (E ■ R) 
=o 



-eoc\Ef- 



(3.51) 



(3.52) 



Pour avoir une expression explicite du vecteur de Poynting, il nous reste a 
evaluer \E\. Pour cela, choisissons un cas simple : 
Pour la particule 

= R 2 (t) = 
R 3 (t) = R(t) 

L'observateur se trouve en 

r x = x\ = r 
y = z = 

Nous appellerons 9 l'angle entre l'axe Oz et la direction du point d' observation. 
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Avec ce choix de parametres, nous avons les relation suivantes 

v x v = ; \R x v\ = Rvs'mO 



d'oii 



R x 



[R+- C R) xv 



R vsin9. Ce qui donne pour \E\ 



R 2 vs'm9 

^0 C 2( R+ ^R 



, ^± t O I L U DILI V , . 



Nous reportons ce resultat dans [3.52] et nous avons : 

Re 2 1 1 . 9. 1 



S = — enc— —r ^^R 4 v 2 sin 2 9- 

i?(4vre ) 2 c 4 R 6 [1 



II cos f 



Soit : 

q = e " f 1 i> 2 sin 2 fl 

16vr 2 e c 3 ^ ii J R 2 h _ ^osfl] 6 1 j 

De la, nous pouvons calculer tP, l'energie emise par le systeme, en integrant le 
vecteur de Poynting sur une surface : 

dE e 2 v 2 [ w sin 2 9 d(cos9) 

—— = b ■ as = ?, 5- / «-27T 

dt J 167T 2 e C 3 Jq \^_vcgse^ 



Ce qui donne 

in 9.-9. 1 i 

5? 



dE e 2 v 2 l+ lv " 



dt Qireoc 6 A 



c 1 



(3.55) 



Ce resultat montre que si v — > c alors ^ — > 00. Ce qui est impossible, nous en 

deduisons done l'impossibilite d'accelerer une particule de t> — > c. 

Prenons maintenant la limite non-relativiste de la formule [3.55], ce qui revient a 



faire - — > 0. Nous obtenons alors : 

c 



dE e 2 v 2 p 2 



dt 67T£oC 3 67T60C 3 

Ce qui est exactement la formule du dipole obtenue en [2.61]. 



(3.56) 
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3.2.6 Structure de la radiation 



Regardons plus en detail la structure de la radiation lorsque v — > c : d'apres la 
formule [3.54], nous pouvons recrire la dependance du module du vecteur de 
Poynting en fonction de Tangle 6 



W)\ 



sin 



'1 



cos 



1 



V I V 

c "r 2c 



2 ) 



6 



Cette approximation pour petit, nous donne l'expression : 

2 



|5(0)| 



Cherchons 6q, la position approximative du maximum : 

df 1 e 2 
— = . - 3 

d0 



'1 



fl2 



1 



1 - - + — = 3<9 2 => = 1 - 

c 2 2 

Ce qui donne un maximum lorsque #q ~ y^f (l — \ ) 



(3.57) 
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Chapitre 4 



Equations de Maxwell et 
transformations de Galilee 



4.1 Principe d'invariance et proprietes de 

l'espace-temps en dynamique Newtonienne 

Considerons un probleme typique en dynamique Newtonienne : la description du 
mouvement de N particules, ayant des interactions de paires, dependant 
seulement des distances entre ces particules. Nous pouvons ecrire l'energie de ce 
systeme 



avec E cin = J2i^ et E pot = J2i,j i^j U (\xi - xj\). L'Hamiltonien (ici l'energie) 
est invariant sous les transformations du type : 
• Translation : 



E — E c i n -\- Epot 



,2 



Xi > Xi -\- Xq 



Nous savons qu'il existe une relation profonde entre symetrie et loi de 
conservation. Dans le cas present, la symetrie mise en jeu implique la 
conservation de Pimpulsion totale P = YliPii avec ^ a definition habituelle de 
l'impulsion pi = rriiXi. 
Nous allons verifier cette conservation : 





E dU dxj __ 
i dxi dxg 




au 

i dxi ' 



• Rotation : 



OijXj , OO t = 1 



3 
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Dans le cas de la rotation, nous savons que la quantite conservee est le moment 
cinetique : 

L = ^2 Xi x pi 

i 

Pour le verifier, nous pouvons etablir le meme type de preuve que plus haut. 
• Translation dans le temps : 

t -► t + t 

Dans ce cas-ci, la quantite conservee est l'energie. 



4.2 Electrodynamique et transformations de Galilee 

II existe aussi des invariances plus subtiles, les transformations de Galilee en sont 
un bon exemple : 

t^t' = t 



4 



La distance entre deux points est conservee : Xi — Xj 
L 'acceleration est conservee : Xi — > x^. 
Du point de vue du sens physique, une telle transformation revient a faire un 
changement de referentiel, ce que nous pouvons representer par le schema 
suivant : 



at=0 
Ox=Ox' 



► x 



Vitesse v 



Nous sommes ammenes aux conclusions suivantes : 

i) La physique est la meme dans tous les systemes de coordonnees en 
mouvement relatif constant les uns par rapport aux autres. 

ii) Le temps est absolu. 

iii) II n'existe pas de mouvement absolu. 

Nous allons maintenant essayer d'appliquer a l'electrodynamique les 
transformations de Galilee. Que se passe-t-il pour la densite de charges p et pour 
la densite de courants j ? 

p est, par definition, la charge par unite de volume, comme l'element de volume 
d 3 x ne change pas, p reste invariant sous la transformation de Galilee. 
La densite de courant, par contre, change : j — ► j' = j — vp. 
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Nous avons vu, en electrostatique et en magnetostatique, que les potentiels 
scalaires et vectoriel sont relies a p et j par les relations : 



1 f p(x')d 3 x' 
47reo J |a? — a?' | 



A 



4tt€qc 2 I \x~x'\ 

les potentiels sont done transformed comme suit 



j(x')d 3 x' 



A^A' 



A 



Ainsi p et 4> se transforment comme le temps, alors que j et A se transforment 
comme une coordonnee. Ce que nous pouvons mettre sous la forme schematique : 



t <-> p <-> 



2 A 



Ayant ceci en tete, nous allons prendre l'equation de Maxwell pour le potentiel 
scalaire, et verifier son invariance : 

<r ot z eo 

Nous prenons dans cet exemple le cas de deux referentiels en translation l'un par 
rapport a l'autre le long de l'axe des x, soit la transformation 



' x' 



vt 



y = y 

z' = z 

t' = t 



II nous reste a evaluer comment les operateurs vont se transformer 



a dx' _j_ a at' 



a_ 
at 



a_ _ _ 

dx dx' dx ~ r at' dx 

_ _a_dx^ , a_d£ _ a 
- dx' at at' at 



a 

dx' 



at' v ax' ' 



ce qui donne : 



a 2 
W 



a 2 



a^_ 

dx 2 



a 2 

" dx' 2 
2 V JLJL _i_ 
zu at' dx' + u 



2 a 2 

dx 71 



Nous en arrivons a l'expression suivante pour l'equation de Maxwell dans le 
nouveau referentiel : 



-Ay 



i ay 



2v d d 
c dt' dx' 



2 fl2i/ 



+ 



V & 



c? dx' 2 



P_ 



Les equations de Maxwell ne sont done manifestement pas identiques suite a une 

transformation de Galilee. 

Trois hypotheses s'offrent a nous : 
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i) II existe un mouvement absolu. 

ii) Les equations de Maxwell sont fausses. 

iii) Les transformations de Galilee ne sont pas les bonnes. 

II s'avere que les hypotheses i) et ii) sont exclues par les faits experimentaux. 



4.3 Transformations de Lorentz 

II nous reste done a envisager l'hypothese (iii) du paragraphe precedent. Nous 
allons done chercher une transformation qui laisse les equations de Maxwell 
invariantes. 

Nous restons ici dans le meme cadre que plus haut : seules les coordonnees x et 
t sont affectees par la transformation. Nous nous limitons done a un espace de 
dimension 2. Nous representons alors la transformation cherchee par une matrice 
2x2 notee A, telle que nous ayons : 

ct' \ . ( ct 
= A 



x 



avec 

A 



An Ai 2 
A21 A22 



Nous obtenons, en terme de coordonnees, les equations : 



ct' = Anct + A12X 
x' = A 2 ict + A22X 



Les operateurs se transforment : 



d a d 1 a 19 

qz. — A22 -foj + A12 - op 

cat ~ Ji21 dx' iin c at' ' 



ce que nous pouvons ecrire sous forme matricielle : 



eft = ( iVn 1V21 ( c|F =A T ( cf 7 

m J V A i2 A 22 J V m J \ m 



Nous remarquons que l'operateur — + ^2^2 qui intervient dans l'equation de 
Maxwell, peut etre ecrit comme : 



Id d \ ( 1 



13 

rJJi Jh-j { -I A |f. 



II se transforme done comme suit 



id a\,/i 



cdtf dx' \ -1 



dx 1 
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Ainsi, si 



Al J _°! J _° ] ,4,, 



alors les equations de Maxwell sont invariantes. L'equation [4.1] peut s'ecrire en 
terme de coordonnees, ce qui donne trois equations independantes : 



A2 _ a2 _ i 
7i n 7i i2 — 1 
A 2 _ A 2 - _i 

yv 21 7V 22 — 1 

A11A21 - A12A22 = 



Nous avons quatre inconnues et trois equations, nous pouvons done introduire un 
parametre pour obtenir une famille de solutions. Appelons x ce parametre, nous 
exprimons alors la solution comme : 

An = coshx A12 = — sinhx 
A 2 i = - sinh x A22 = cosh x , 

oil x es t appele la rapidite. Nous avons pour les coordonnees : 

ct' = ct cosh x ~ x sinh X 
x' = —ct sinh x + x cosh x 

Quelle est ^interpretation physique de sinh x ? Prenons x' = ceci nous donne 
* = ^ = t*. Alors, 

v 

tanh x = ~ = 
c 

cosh x = — 1^=^= = 7 sinh % = — = 7/? 



et la transformation a la forme suivante 



x' = j(x - (3ct) ; ct' = j{ct - fix) . (4.2) 

Pour obtenir la transformation inverse, il suffit de faire (3 — > —/3 

Conclusion : la physique ne change pas si nous changeons de referentiel, et si 

nous changeons les coordonnees suivant la formule [4.2]. 

De telles transformations sont appelees des transformations de Lorentz. Nous 
pouvons generaliser [4.2] pour une direction arbitraire de la vitesse v : 

xj| = 7(xj| — (3ct) ; x' ± = x± ; ct' = y(ct — (3 ■ x) (4-3) 
oil x\\ = j£ (3 et x± = x — xu . 
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Les transformations correspondantes pour p, j, <j>, A sont donnees par : 



cp 



CA '\\ 

A> 



l[cp- P-j 
7 ~ Pep) 
3± 

1 U-p- (cA) 

7 (ciy - p<p) 

A± ■ 



(4.4) 



On conclut done que (cp,j) et ((f), cA) transforment comme (ct,x) sous les 
transformations de Lorentz. 

Nous voyons aussi que les transformations ne sont valides que pour P < 1. 



4.4 Potentiel d'une charge ponctuelle en mouvement 

Pour verifier notre travail, nous allons calculer le potentiel d'une particule 
chargee se deplagant a une vitesse constante v dans le referentiel du laboratoire. 
Partons de : 

6 = 6 ± 

V ^ Ane r U 5) 

,4 = V ' ' 

Ainsi, a partir des transformations de Lorentz nous avons : </)' = 7^ soit 
j.i e 7 1 



47re 



Une autre possibility s'offre a nous pour faire ce type de calcul : utiliser les 
potentiels de Lienard-Wiechert : 

Pour cela, nous allons nous placer dans le systeme de coordonnees pour lequel la 
charge electrique se deplace a une vitesse — v : 



47re |£_£| + ^^) 



c 



ou le temps retarde t est donne par 



t' = t -\ E —^ \R-x\ = c(t-t') . 

c 



Or, comme R = —vt', nous obtenons 

t' = t--^/y 2 + z 2 + (-vt' -x) 2 
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done : 



\R-x\ + — y — = c(t - t') + -(vt' + x) = c 

Nous devons trouver t' : 



vx. ., v 2 J 

(t+^)-t'(l-- 2 ) 



(y 2 + z 2 + (vt' + x) 2 ) = c 2 (t - t') 



l\2 



De la, nous tirons 



l-^)t' = t 



vx 1 



done : 



v-(R-x) 



(x + ^) 2 + 1 - % )(y 2 + z 2 ) . 



Ceci n'est rien d'autre que le denominateur du potentiel cj), en le reportant nous 
obtenons : 

e7 1 



cj> = 



47re ° ^y 2 + z 2 + ^(x + vt) 2 ' 



ce qui est exactement l'expression obtenue a partir des transformations de 
Lorentz. 

Pour le potentiel vecteur, nous avons : 
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Chapitre 5 



L 'elect rodynamique avec les 
quadrivecteurs 



Nous arrivons maintenant a une contradiction. D'apres les equations de Maxwell, 
nous trouvons les transformations de Lorentz. Mais, en partant des equations de 
Newton, nous obtenons les transformations de Galilee. 

Nous en deduisons que l'equation de Newton doit etre fausse pour de grandes 
vitesses v. 

Notre but est done de trouver une generalisation des equations de Newton, qui 
soit coherente avec les transformations de Lorentz. Pour y parvenir, nous allons 
developper un formalisme 4-dimensionnel pour les equations de Maxwell. 

5.1 Definitions : quadrivecteurs covariants et 
contravariants 

Definition : un vecteur, dans l'espace a 3 dimensions, est un ensemble de 3 
nombres qui se transforment comme (x, y, z) sous une rotation du systeme de 
coordonnccs. 



Les operations suivantes determinent des grandeurs vectorielles ou scalaires : 



A = A = A 2 ; A' = OA ; 00 T = 1 




nous avons des invariants : le produit scalaire 

A ■ B = A T B -► A T T OB = A T B 



< 



V0 vecteur 
V • A scalaire 
Ax B vecteur 
VxB vecteur 



v 
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IMPORTANT : Pour la suite, nous utiliserons les notations suivantes 
x° = ct ; x 1 = x ; x 2 = y ; x 3 = z, soit sous forme vectorielle : 

( ct\ 

X 

V 

V * J 

ou bien x M , avec [i = 0, 1, 2, 3, en mettant bien l'indice /U en haut (quadrivecteur 
contravariant, voir ci-dessous). 

Par la suite, nous ne considererons que les transformations de Lorentz avec 
v || e* x , c'est-a-dire les 'boosts' suivant x. 

De cette fagon, nous pourrons ecrire les transformations de Lorentz sous une 
forme elegante : x — > x' = Ax, avec la matrice A definie par : 



A 



/ 7 -7/3 \ 

-7/3 7 

10 

\ o o 01/ 



(5.1) 



Definition : un quadrivecteur contravariant est un ensemble de quatre nombres 
qui se transfer me nt comme x^ si nous changeons de referentiel. 
Exemples : = (cp,j) T , = (<j),cA) T . 

Remarque : une transformation peut etre ecrite en terme de composantes 



comme : 



La matrice A satisfait la relation A T GA = G, ou la matrice G est : 
G = diag(l, -1,-1,-1). 

Definition : La matrice G est appelee la metrique ; ses composantes sont g^ u , 
avec les indices en bas. 

Definition : Un quadrivecteur covariant est un ensemble de 4 nombres qui sont 
derives d'un quadrivecteur contravariant a M = (ao,a) T comme suit : 



a 



(a , -a) 



Par la suite, nous utiliserons la convention d'Einstein 

3 



V _ \ " V 

9^a = 2^ 9^a , 



u=0 



c'est-a-dire le symbole pour la somme est supprime avec la convention qu'on doit 
toujours sommer sur des indices repetes en haut et en bas. 
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Notons ici la transformation d'un quadrivecteur covariant : 
Sous la forme matricielle l'equation precedente est ecrite 

(i)- a (S) - ga 00 - gao {-^ <^ r (. 3 )- 

Ainsi, nous pouvons calculer le 'produit scalaire' 1 comme : 

a!^ = A/a a A^ = g%ay = SfiaJbP = a a b a (5.2) 

ou 

(«o , • ($) = («o , -a) A-A = (ao , -a) • ( 6 |) . (5.3) 

Les produits scalaires de quadrivecteurs ne changent pas si nous changeons de 
referentiel, done aussi la norme 

x » XfM = ( x ) 2 - \x\ 2 = c 2 t 2 - \x\ 2 (5.4) 

est un invariant. 

5.2 Les equations de Maxwell en 4 dimensions 

Maintenant que nous avons decrit le formalisme lie au quadrivecteurs, nous 
allons pouvoir exprimer l'electrodynamique, et en particulier les equations de 
Maxwell, sous une forme adaptee a l'etude de problemes relativistes. 

5.2.1 Exemple 1 : 

Regardons queries sont les proprietes de transformation de 

d 
dx^ 

Nous avons deja trouve que : 

( ¥ ) = A T ( W ) 5 on (4°) = (A-T (i) • (5.5) 

V 83 J \ dx> J \ dx> ) \dx J 

On a mis des guillcmets, parce que ce n'est pas un produit scalaire d'apres la definition 
mathematique. La raison e'est que la 'metrique' g M „ n'est pas definic-positive, mais indefinie, 
e'est-a-dire il y a des quadrivecteurs non-nuls avec la norme nulle et meme des quadrivecteurs 
avec la norme negative. 
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Done = est un vecteur covariant. Ainsi pour le quadripotentiel A 11 , nous 
pouvons calculer le scalaire suivant : 

ce qui est exactement la forme de la condition de la jauge de Lorenz que nous 
avions pour les equations de Maxwell. Nous en deduisons que la condition de 
Lorenz ne change pas si on change de referentiel. 



5.2.2 Exemple 2 : 

Maintenant, regardons Poperateur de d'Alembert : 

D = ()„()'' - I iH-A 



c 2 dt 2 

Le produit scalaire de quadrivecteurs est invariant, or les equations de Maxwell 
sont de la forme : 

UA»=— f (5.6) 
ceo 

La conservation de la charge electrique : 

d^f = (5.7) 

est un invariant de Lorentz. 



5.2.3 Transformation des champs, tenseur champ : 

Nous devons comprendre comment les champs electrique et magnetique se 
transforment. Pour cela, nous allons considerer la quantite suivante : 

= d^A u - dyA^ (5.8) 

Ceci est appele le tenseur champ ; il est labele par deux indices. Nous donnons ici 
quelques proprietes : 

i) F^y = —F UfJi , le tenseur champ est antisymetrique. 

ii) Le tenseur champ est invariant sous les transformations de jauge 

A^ ► A^ = A^ — d^a . 

Preuve : 

F^v -> d^(A u - d v a) - d^A^ - d^a) = _ Q , 

Fftv - dpd v a + dyd^a = F ilv 



42 



iii) Les composantes du tenseur champ sont reliees aux champs electriques et 
magnetiques : 



plV 



I 

V 





-E x 

~Ey 

-E K 



E x 

cB z 
-cB„ 



Ey 

-cB z 




E 

cB, 



\ 



v 

cB, 



/ 



exemple de calcul des composantes typiques : 

1 



Foi = d A l - diA 
F 23 = d 2 A 3 - d 3 A 2 



d t {-cA x ) - d x (f> = E x 
d y (-cA z ) - d z (-cA y ) -- 



-cB r 



(5.10) 



Grace a ce tenseur, les equations de Maxwell peuvent s'ecrire sous cette forme : 
Equations inhomogenes : 



= — f 

ce 



V-E = f 

c 2 VxB =^ + # ' 
eo at 



(5.11) 



ou 



= d ll A u - d u A^ 





—E x 


— Ey 


-E z \ 


E x 





-cB z 


CBy 


Ey 


cB z 





-cB x 


\Ez 


-CBy 


cB x 


o / 



Equations homogenes : 



c <~Jv L pa 



V-B = 

VxE+f =0 



(5.12) 



ou e^f" 7 est le tenseur totalement antisymetrique (tenseur de Levi-Civita) . 



uvpa 



1 si (fivpa) une permutation paire de (0123) 

— 1 si (fivpa) une permutation impaire de (0123) . 

sinon 



De la, nous voyons une fois de plus que les equations de Maxwell ne changent pas. 
La conservation du courant peut etre exprimee comme suit : 



d v d p 



symetriquc 



antisymetrique 



1 

ce 



duf 







(5.13) 



_ g m g v0 Fa/3 
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Transformations de Lorentz de E et B : 

Pour trouver les transformations de Lorentz de E et B, il est maintenant 
sumsant de trouver la transformation de F^ u . 

F 1 ^ = A^ a A u /3 F al3 

prenons, par exemple, la vitesse parallele a l'axe (Ox), alors, nous avons : 



A = 



( 7 -h o \ 

-07 7 

10 

V o 001/ 



Soit, done : A° = 7, A° 1 = -£7, = -^7, A\ = 7. 
Nous calculons les transformers : 

F' 01 = A aA 1 /3 F al3 

= A^A^F 00 + A° 1 A 1 F 10 + A° A\F 01 + A^A 1 ^ 11 

= p^ F W + 7 2 F 01 = F 01 7 2 (1 _ 02 ) = F 01 

Nous en deduisons done, pour le champ electrique : 

^ E' x = E x 

Nous faisons un second calcul pour la composante perpendiculaire : 



F' 02 = A° a A 2 pF al3 
,0 7702 



77«2 



A U ^F 



7712 



d'ou 



A U F + A\F 



-E' y = j(-E y + (3cB z ) . 



(5.14) 



(5.15) 



(5.16) 



Nous trouvons une equation analogue pour la troisieme composante. Pour 
champ electrique, nous pouvons done ecrire : 

#11 = E\\ 

E' ± = 7 (E ± -cB xf3) . 



(5.17) 

(5.18) 

(5.19) 
le 

(5.20) 
(5.21) 



Nous effectuons le meme type de calculs pour le champ magnetique : 
F' 23 = A 2 a A 3 pF a/3 = F 23 ^B' X = B X 



F' 12 = A 1 a A 2 /3 F al3 = A\F a2 

= A\F 12 + A^F 02 = 7F 12 - 7/3F' 



02 



(5.22) 
(5.23) 
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ce qui donne : —cB' z = j(—cB z + f3E y ). Done : 



B'\\ = B\\ (5.24) 
B'± = l{B± + X -E x $) . (5.25) 
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Chapitre 6 



Les particules dans les champs 
electromagnetiques 

6.1 Les equations de Newton : Cas relativiste 

Pour discuter du mouvement des particules, nous devons ajouter aux equations 
de Maxwell les equations de Newton, dans lesquelles la force de Lorentz va 
intervenir. Pour un cas non-relativiste elles s'ecrivent : 

m 3 -^- = q [E + v x B) (6.1) 

Nous voulons, dans cette partie, trouver une generalisation des equations [6.1] 
pour la mecanique relativiste. L'idee est la suivante : nous allons ecrire les 
parties gauche et droite de [6.1] sous forme de quadri-vecteurs. 

1. Considerons, dans un premier temps la partie gauche de Pequation [6.1] : 
La vitesse 4? a des proprietes de transformations compliquees. Nous allons 
done essayer de trouver une quantite plus simple. 

est un quadri-vecteur dx^ est done egalement un quadri-vecteur. La 
quantite dx^dx^, appelee l'interval, est un scalaire. Nous notons : 

ds 2 = dXf.dx' 1 = g^dx^dx" (6.2) 

Si nous considerons une particule en mouvement, avec les coordonnees 

x* = x\t) , i = 1,2,3 => dx i = v l dt 

Nous avons alors : 

ds 2 = c 2 dt 2 - v 2 dt 2 ds = dt\Jc 2 - v 2 = c— 

7 

ds est un scalaire, dx^ 1 un quadri-vecteur, done est un quadri-vecteur. 
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La quantite : 

ds 

est appelee la quadri-velocite (ou quadri-vitesse) , dont les composantes 
sont : 

c— — = = j(c, v) 
ds 

avec v^v^ = 7 2 (c 2 — v 2 ) = c 2 . Ainsi, pour avoir un quadri-vecteur dans la 
partie gauche de [6.1], nous pouvons ecrire : 

2 d 2 x' 1 . , 

mc ^ (6.3) 

qui est si v -C c. 

2. Considerons maintenant la partie droite de [6.1] : 

cette partie est composee de : E + v x B. Essayons une combinaison du 
type : 

F^ v v v pour pi = 1,2,3 
=> F 10 uo + F 12 v 2 + F 13 v 3 = -F w v° - F 12 v 2 - F 13 v 3 
= "fEic + jcB 3 v 2 — jcB 2 v 3 = 7(.Ei + (v x -B)i)c 

Done : 



F 4 % = 7 + (v x B) 1 J c (6.4) 

Ceci est la force de Lorentz a un facteur 7c pres. 

En mettant les equations [6.3] et [6.4] ensemble, nous obtenons la forme 
relativiste de [6.1] qui s'ecrit : 

mc 2 ^ = -F^Vv (6.5) 
ds z c 

Cette derniere equation est automatiquement invariante sous transformations de 
Lorentz. Elle est identique aux equations non-relativistes dans la limite des 
petites vitesses v <C c. Notons aussi que nous obtenons une equation 
supplementaire pour fj, = 0. 
Ecrivons chaque composante explicitement : 

1. Composante temporelle : 

mc 2 d ( c \ d q 0i 

' ' — 7— [m'jcj = -t 1 Vi 



yjc 2 - v 2 dt V Vc 2 - v 2 J dt 

F 0i Vi = F 0i v l = 1 {E-v 
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ce qui peut se mettre sous la forme : 

d I mc 2 



, , , = qE ■ v 

dt \ f[Z^ J *-v-' 

, V <? / travail 



L'energie est donnee par £ = , mc2 = 7717c 2 , done ^§ = o-E • v. Notons que 



pour f = 0, nous trouvons 



£ = mc 2 

qui est la formule d 'Einstein. 
2. Pour les composantes spatiales : 



mc 2 7-^- I — = I = - jc(E + vxB) 



cdt I _ A ~^fl ] c 

soit : 





V 2 J 

C 2 / 


dp 
dt ~ 


q(E + 


p = 


mv 



Ceci nous permet d'introduire le quadri-vecteur impulsion : 

= m?j M , p° = —, p % — »■ p 
c 

avec 



£ = y^m 2 r?~+~?^c 2 
jj = 7777-1/ 

6.2 Applications pour des cas simples : 

6.2.1 Champ electrique constant : 

Prenons le cas d'un champ electrique oriente suivant l'axe (Ox). 

f E x ^0 
\ E y = E z = 

Prenons egalement vf n = 0. 

Nous avons alors les equations suivantes : 

J ^ = 1 E 

) _ dp Z _ p, 

l dt ~ dt ~ " 



(6.6) 



(6.7) 
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dont les solutions sont 



p x = qEt + p\ 
Py=Py 



(6.8) 



P expression de Penergie est done : 



£ = ^m 2 c 4 + ( P 0)2 C 2 + ( q Et + p° x ) 2 c 2 

pour la vitesse : 



2 p [(qEt + p Q x )e x + p y e y ] 

v = c — — 



c 2 



/ m 2 c 4 + ( p 0)2 c 2 + ( qEt + p 0)2 c 2 

Notons que pour c — > oo, nous trouvons : 

^ 2 p [(9^+^)4+^] c 2 / g£A _ ^ ^ gEt 

w = c z - = = \v u x H e x + v% = v H 

£ mr \ m / y m 

qui est Pexpression de la vitesse non-relativiste. 

Pour t — > oo : u — > c. Nous trouvons pour la trajectoire : 

dt ^/m2 C 4 + (pO)2 C 2 + (g^+pO)2 c 2 

rfy =r 2 l °' yJ 

<** ^/m 2 c 4 + (p° ) 2 C 2 + (gfit +p ) 2 c 2 

Ce qui donne comme solution pour x (avec une certaine condition initiale) : 

x = ^^m 2 c 4 + (p° y ) 2 c 2 + (qEt + p0) 2 c 2 
Done si i — > oo, x = ct 

6.2.2 Champ magnetique constant suivant (Oz) : 
Dans ce cas, nous avons Pexpression : 



dp 
~dt 



= q(v x B) 



— £ = qE ■ v = f = const , 7 = — = const 

* ^1 



— = w x w B , = = — — 

at 7m c 
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Nous obtenons les equations suivantes pour la vitesse : 

dv z — n 
dt ~ u 

it = V v"B 
dv v 



Qui ont pour solution : 



Vx 



(6.10) 



(6.11) 



Ce qui donne pour les coordonnees (avec une certaine condition initiale) : 



z = v z t 
~cos{u B t + 



> y = Z^ sin(w B t + <j)) 



Ce qui donne une spirale comme suit : 




(6.12) 



Pour v z = 0, nous avons le rayon qui est : rs 



V 



6.2.3 Formule de Larmor relativiste 

Nous connaissons deja une fagon d'obtenir la radiation d'une particule acceleree : 

1. utiliser les potentiels de Lienard-Wiechert pour E et B, 

2. calender le vecteur de Poynting a une distance R assez grande, 

3. integrer le vecteur de Poynting par rapport aux angles et au temps pour 
trouver l'energie totale emise par la particule. 

Maintenant, nous pouvons, grace a la mecanique relativiste, trouver une 
methode beaucoup plus simple. 
Donnons ici l'expression classique : 



7 2 -u 2 



d£ = ^dt 



(6.13) 



l'energie est la composante du quadri-vecteur 
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le temps est la composante du quadri-vecteur x^. 

Or est un invariant, done une generalisation de [6.13] doit etre invariante de 
Lorentz. 

Si nous remplacons v par ^pr, nous avons : 

d£ q 2 ( d 2 x^ 2 d 2 x^ 2 \ 

" c ~ZT7T C ( 6 - 14 ) 



dt 67reoc 3 V ds 2 ds 2 



scalaire 

Nous voyons que la limite non-relativiste reste correcte. D'autre part, nous 

2 d 2 x^ q 



pouvons remplacer ^-Mr par dans [6.14] grace aux equations du mouvement : 



ds z mc 



Nous obtenons alors : 

d£ 1 



F^v u F w vP (6.15) 



dt Gireoc m 2 c 4 v MP 

6.3 Les accelerateurs 

II s'agit ici de reprendre les calculs effectues a la section ci-dessus, et de les 
appliquer aux cas des accelerateurs. 

6.3.1 Accelerateur lineaire 

Dans le cas d'un accelerateur lineaire, nous avons E \\ v. Nous pouvons done 
ecrire : 

~ = [F 01 v 1 F 01 v 1 + F w v F w v°] = -t^L^E 2 (v 2 - c 2 ) 7 2 

en simplifiant, nous obtenons : 

d£ q A E 2 
dt 67reom 2 c 3 

qui ne depend pas de Penergie de la particule ! 

6.3.2 Radiation synchrotron 

Pour un synchrotron, nous sommes dans le cas oil le champ magnetique B est 
perpendiculaire a la vitesse de deplacement de la particule. Dans ce cas, le calcul 
est le suivant : 

1 rz7>12 1-1 2 , t-i21 j-i 11 9 4 2 D 2 2 

-tt = ~a [F v 2 F 12 v + F AL v 1 F 21 v L \ = K-^rfB v 
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Ce qui se simplifie en : 



d£ 
~dt 




67reom 2 c 3 ( 1 




Dans cette derniere expression, le facteur souligne nous rappelle que la radiation 
augmente lorsque v — > c. Ce qui implique que les accelerateurs circulaires ne sont 
pas efficaces pour de hautes energies. En effet, la particule pert de Penergie en 
emettant des photons dans une direction tangente a sa trajectoire. Notons 
toutefois que la radiation synchrotron a des caracteristiques particulieres qui la 
rendent interessante pour d'autres applications. 
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Chapitre 7 



Champs electrique et 
magnetique dans un milieu 



7.1 Les equations de Maxwell macroscopiques 

Notre objectif, dans ce chapitre, est de deriver les equations de 
l'electrodynamique macroscopique (qui ont ete donnees par Lorentz en 1902). 
Pour ce faire, nous partons des equations microscopiques : 



ou nous avons les champs microscopiques qui sont represented par des lettres 
minuscules. Nous rappelons egalement la relation : = c 2 . Pour les milieux 

continus, e et b sont, en general, des fonctions compliquees de x et de t a des 
echelles tres petites : 

La distance typique de variation est la taille de l'atome ~ 10 _8 cm, le temps 
typique est done de l'ordre de 3 ^oiocm/s ~ lO -18 *- 

Experimentalement, les ordres de grandeurs sont : t » 10 _18 s ; L » 10 _8 cm. 
II nous faut done introduire des quantites macroscopiques : 



ou f{x) est une fonction qui tend vers zero lorsque |x| > R, avec R grand devant 
les distances inter- atomiques, mais petit devant la taille de la region 
d'observation. 

Nous donnons ici des exemples pour la fonction / : 



macroscopique 




(7.1) 



E{z) = J d 3 x e(x)f(x- z) = {e(z}) 
B(z) = fd 3 x b(x)f(x-z) = (b(z)} , 



(7.2) 
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1. 

2. 
3. 



4^ > r<R 

, r> R 



(ttR 



2^-3/2 e _ r 2 /i? 2 



d 3 x f(x) = 1 



(7.3) 

(7.4) 
(7.5) 



si R — > 0, nous avons une fonction J. 
Donnons ici quelques relations importantes : 
En utilisant l'integration par partie 

Vg{e{z)) = fd 3 xe(x)V?f(x-z) = 
J d 3 xe(x) (-Vsf{x- z)) = J d 3 xV 3 ef{x- z) 



(V-e) 



alors 



D' autre part, nous avons 



V-(e): 

m 

dt 



(V-e") . 



K dt' 



Pour obtenir les equations de Maxwell macroscopiques, nous devons alors 
une moyenne des equations [7.1] : 



(7.6) 

(7.7) 

(7.8) 
faire 



V-b 

V x e = - 
V • e = 



dt 



= 
-2- = 



> V-B-- 
■ V x E 
V-E = 







9B 

at 



(7.9) 



-^-V x 5 

£0/^0 



ag 
at 



01 



II nous reste maintenant a determiner ce que sont les quantites (rj) et (j), 
c'est-a-dire les moyennes de la densite de charges microscopique r\ et de la 
densite de courant microscopique j. 



(i])=p{x,t)-V-P 



(7.10) 



ou nous avons note : p la densite de charges macroscopique (densite de charges 
libres) et J la densite de courant macroscopique (densite de courant de charges 
libres). Par exemple, ce sont des electrons ou des ions qui peuvent se mouvoir 
librement. P est la polarisation macroscopique, qui peut etre definie comme la 
densite de moment dipolaire p = ed, done 



V ^ 



Vr, 
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De la meme facon, M est l'aimantation macroscopique, donnee par 



n 

avec rh n = I n Sn = |a? n x v n . 

Demontrons que rh est exactement donnee par la relation standard rh = ISn, oil 
S est l'aire, I le courant total et n est le vecteur normal a la boucle de courant. 
Le courant total est donne par : 

Aq 
At 

Pour un mouvement circulaire, definissons combien de fois nous croisons un 
point donne sur la boucle pendant le temps At. La distance est vAt, done le 
nombre de fois est 

vAt 



2ttv 



la charge totale est done : 

qvAt qv 

=4> done le courant 

2irr zirr 

La norme du moment magnetique est alors donnee par : 



tq 1 V 2 1 

lb = irr = -qvr (7.11 ) 

2irr 2 



ce qui correspond exactement a : 



rh = —x x v (7.12) 
2 v ' 



7.2 Champs H et D 

Determinons maintenant le sens physique des differents termes : 

1. Terme V • P : 

Considerons des molecules dipolaires du meme type que celle de la partie 
gauche de la figure [7.1]. Si ces molecules sont distributes de manieres 
chaotique, il n'y a pas de charge induite, mais si elles sont ordonnees, 
comme dans la partie droite de la figure, la densite de charges induite 
Ppoi = —V • P donne une fonction non nulle du type de celle de la 
figure [7.2] 

2. Terme ^ : 

Reprenons la meme distribution que celle de la figure [7.1], mais ajoutons 
un processus dynamique : une inversion de la polarisation des molecules 
(par exemple une rotation), les charges + vont en bas et les charges - vont 
en haut. Alors, il est clair que le terme ~ j est une densite de courant. 
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Vide 



r. t 
i: i: 



Vide 

Fig. 7.1 - Molecule dipolaire, et une distribution possible 

P(z) <ip 



Fig. 7.2 - Distribution de P et (jj) pour la fig. 7.1 



3. Terme V x M : 

Ce terme represente une densite de courant de surface. 

Pour tenir compte de ces termes physiques que nous venons de voir, 

nous introduisons ici des nouveaux champs : Le champ magnetique H 

et l'induction electrique D. 1 Les definitions pour D et H sont donnees par les 

equations 

f D = e E + P 



H 



Ho 



-B - M 



(7.13) 



(7.14) 



Ces nouvelles definitions impliquent une reecriture des equations de Maxwell 
comme suit : 

( Vx#-| = J, V-D = p 
( Vx£ + | = 0, V-B = 

7.3 Solution des equations de Maxwell 

Comment resoudre les equations [7.14,7.13] ? En effet, nous devons determiner 
les champs E, B, D, H, nous avons done deux inconnues de plus que pour les 
equations dans le vide. Les solutions ne peuvent etre trouvees que s'il existe des 
relations entre elles. 

Les relations les plus simples possibles sont pour des milieux continus isotropes, 
nous trouvons alors : 

\ D = eE . . 

) \b (7 - 15) 



H 



1 Pour eviter une confusion entre H et B on appelle le champ B l'induction magnetique. 
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avec la permittivite e et la permeabilite \i qui sont des constantes 
caracteristiques a chaque materiel, e et \i sont determinees soit 
experiment alement, soit grace a des calculs de microphysique. 
Les equations [7.15] se generalisent pour un milieu anisotrope : 



A 

Bi : 



- €ijEj 



(7.16) 



ou eij et fiij sont des tenseurs. 
Remarque : 

les relations lineaires les plus generales possibles sont donnees par : 

Di(x,t) = J d 3 x'dt' ^(x^x* ,tf)Ej{x* J) 
Bi(x,t) = fd 3 x'dt' iM j (x,t,x',t')H j (x*,t!) 

D 'autre part, pour des champs intenses, des termes non-lineaires peuvent 
apparaitre. 



(7.17) 



Fig. 7.3 - distribution quelconque de molecules dipolaires 



Dielectricite : 

Construisons maintenant la theorie microscopique d'un dielectrique. Nous 
connaissons l'energie du dipole en champ electrique : 

£ = —E ■ p = —pE cos 9 



ainsi, certains dipoles sont energetiquement favorises, comme c'est le cas dans la 
figure 7.4. Nous pouvons definir une quantite globale, la polarisation moyenne, 



Fig. 7.4 - Le dipole de gauche est energetiquement plus favorable. 
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Ep cos 6 

en utilisant la probabilite de Boltzmann V ~ e \ 

Ep 

J p e^T cos Od cos 6 d(f> 1 2 

J ekr d cos9 dcf> po jF petit 6 fci 



Nous en deduisons alors : 



D = '" E + -/Tf 



Ce qui nous donne pour la permittivite : 

/ 1 p 2 n 

Paramagnetisme : 

Nous faisons le meme type de calculs pour un moment dipolaire magnetique fa 
et nous obtenons 

1 nm 2 B 

H=^B- ln f B ^ (7.19) 

( nmVoV 1 - ft , nm Vo\ , 7on s 

»=^\ l -^r) ™^\ l+ ^r)- (7 - 20) 

7.4 Les ondes electromagnetiques dans les milieux 
continus 

Nous consider ons ici le cas isotrope, nous avons done : 

D = eE; B = /j,H . (7.21) 

On recrit la loi de Faraday 

VxE+^ = e> (7.22) 
d 2 B dE 

1SF+ Vx-=0, (7.23) 

puis on utilise 

- 3D 1 - dE , . 

Vxif- — = -VxB- £ - = 0, 7.24 
at jj, at 
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ce qu'on introduit en [7.23] et on obtient 



d 2 B dE 



d 2 B 1 / 
§ ^(V(V-VA^) 



(7.25) 



at 2 



i 

C/i 



-AS = . 



ou on a utilise V • B = 0. 

Cette derniere equation est une equation de propagation. La vitesse de la lumiere 
est donnee par : 



1 



v 



1 eoMo 2 e 0^0 , 2 
= c < c 



e/i e //o e/x 



e/i 



(7.26) 



La vitesse de propagation d'une onde dans un milieu est done toujours inferieure 
a c. 

7.4.1 Reflexion d'une onde 

Nous considerons ici le probleme de la reflexion d'une onde electromagnetique a 
l'interface de deux materiaux (cf. figure 7.5). Pour comprendre la physique de ce 





K' 


M-i e i 


/ i 


r' ^ 




V K 






K7 



incident 



reflechi 



Fig. 7.5 - propagation d'une onde a une interface 



phenomene, nous devons tout d'abord trouver les conditions aux bords. Si 
l'interface est tres bien marquee, que devons-nous faire avec les derivees de e et 
de fi ? 
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L'idee est alors d'ecrire les equations de Maxwell macroscopiques sous forme 
integrale, et de les appliquer a l'interface. 

' (1) §DdS = JpdV 



< 




(7.27) 



k (4) jHdl = JJdS+l t jDdS 

II ne nous reste plus qu'a utiliser ces dernieres equations pour un volume et un 
contour specifiques du type de la figure 7.6. 

_^T^ 

(a) 

Fig. 7.6 - (a) petit volume a l'interface; (b) contour autour de l'interface, le 
vecteur t venant vers nous. 

Nous supposons que sur la surface, il peut y avoir des densites de charge 
macroscopiques a (par unite de surface) et de courant macroscopique K (par 
unite de longueur) non-nulles. 

1. Appliquons l'equation (1) de [7.27] a ce petit volume : 

j> D dS = (D 2 - Di)nAS = aAS 

{D 2 - Dt) ± = a (7.28) 

ce qui signifie que la composante normale de D a un saut egal a la densite 
de surface de charge. 

2. Appliquons l'equation (2) de [7.27], nous obtenons le meme type 
d'equation : 

(B 2 - Bi)± = (7.29) 

qui signifie que la composante normale du champ magnetique est continue 
a l'interface. 

3. Appliquons l'equation (3) de [7.27] au contour de la figure[7.6] : 

J Edl = => {E 2 - Ex)\\ = (7.30) 
La composante tangentielle de E est continue. 




02 



4. Appliquons l'equation (4) de [7.27] au contour : 

y (V x H)-dS = j H ■ dl = (fx n) ■ (H 2 - H ± )Al 



J-dS= (K-t) Al , 



ou t est un vecteur tangent a la surface et K est la densite de courant de 
surface. 

Finalement, on obtient une condition au bord pour la composante tangente 
du champ magnetique : 

{H 2 -Hi)\\ = K. (7.32) 

Le saut du champ magnetique H est egal a la densite de courant de surface. 

Nous allons maintenant considerer le probleme de la reflexion dans un cas 
simple : 

<7 = K = 

Nous notons alors l'expression d'une onde sous la forme d'une onde plane : 

E = Eq exp[ik ■ x — iut] (7.33) 

ou la quantite physique est la partie reelle -Ephys = $l(E). 
k est le vecteur d'onde, 
to est la frequence, 

v 2ttv 271 

v u k 

j£ est la vites 
D 'autre part, de la loi de Faraday : 

^+Vx£ = 0, (7.34) 
nous deduisons la forme explicite du champ magnetique : 

B = Bq exp (^ik ■ x — iut 



A = I = ^ = ^ est la longueur d'onde, 
v = j est la vitesse de l'onde. 



-iooBo + ik x Eq = =^ B = —k x E = \k x E 

UJ VK 



(7.35) 



Pour la suite, nous noterons : 

i) onde incidente : Eq, k, u>, v, 

ii) onde transmise : E' , k', to' , v' , 
Hi) onde reflechie : Eq, k", to" , v". 
Nous savons egalement que nous avons 



63 



1. la continuite dans le temps : to = uj' = to" i.e. la frequence ne change pas. 

2. si le plan z = est le plan de reflexion, alors la continuite a z = nous 



donne : 



k ■ x ) = ( k' ■ x ] = ( k" ■ x 

z=0 \ / 2=0 \ J z=0 



Done, les trois vecteurs k, k' , k" se trouvent dans le meme plan. Prenons alors 
k y = 0, ceci implique 

k sini = k sinr = k" smr' 
Comme u = u" et v = v", nous avons k = k" et done : 

sin i = sin r (7.36) 

En plus, nous en deduisons : 



sinr k cuv' v Ui^i n i 



sini k' ujv v y £2^2 ™2 
oil nous avons introduit les indices de refraction : 



(7.37) 



m 



Nous deduisons de l'equation [7.36] que Tangle de reflexion est egal a Tangle de 

Tonde incidente {i = r'), cf. figure 7.7. 

L'equation suivante, tiree de [7.37], est la loi de Snell : 

sini 112 n\ 

= — sm r = — sm 1 . 

sinr ni n 2 

Remarquons que cette equation n'a pas de solution pour ni > ri2 si 

'n 2 " 



1 > iq = arcsm 

\ni 

Nous disons alors qu'il y a reflexion totale pour i > iq. C'est le cas pour une 
interface eau-air avec un angle d'incidence de 43°. 

Pour trouver les relations [7.37], nous n'avons jamais utilise la forme explicite 
des conditions aux bords. Ces dernieres permettent de determiner les amplitudes 
des ondes transmises et reflechies. Considerons un probleme general, le nombre 
total d'equations est 1 + 2 + 2 + 1 = 6 et nous avons 6 inconnues. 
Conditions : 



1. 


(D 2 


-Di) ± 


= 


2. 




-A) 


= 

L 


3. 


(e 2 




= 



exiEo + E'^-ezE'o) - n = 
kx E + k" x E'l - k' x E' ) ■ n = 
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Fig. 7.7 - Schema de principe pour la reflexion-transmission d'une onde 



4. (#2-#i) =0 



1 L(kxE + k»xE»)- 1 Lk>xE> 



x n = 



Considerons pour simplifier, la polarisation lineaire de la figure 7.7 
Solution : 

Supposons que E est parallele a la surface (polarisation TE). C'est-a-dire qu'il 
reste deux inconnues a la place de 6. Nous verifierons cette supposition dans les 
equations par la suite : 

1. Satisfait automatiquement. 

2. Nous avons : 

|&||i£o| sm * + l&'H-E'ol srn ^ ~~ l^'ll^ol smr = 
D' autre part, comme \k\ = \k"\, et |fe|sini = |A/|sinr Nous obtenons : 

| So I + \K\ " \K\ = (7.38) 

3. De meme que le point precedent. 

4. Ici : 

(kEo — kE'o) cos i k'E' cos r 



= 

Mi ^2 
— — (Eq — E n ) cos i cos r = 



Nous utilisons que h = \ jriri d'oii 



fc' V /*2£2 ' 

. I^(E - E%) cos i — \ h^E'o cos r = (7.39) 
V Ml V M2 

Nous avons done 2 equations (7.38, 7.39) pour 2 inconnues (2£q et Eq). 
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Les solutions sont 



E' Q 2ni cos i 

En 



ni cos i + ^ -y/ — nf sin z i 



£// ni cos i — ^ \J n \ ~ n i srn2 i 

En 



n\ cos i + ^J n 2 ~ n i s ^ n 



2 i 



et sont appelees relations de Fresnel. 

Dans la plupart des cas (par exemple pour des frequences optiques), nous 
pouvons considerer que ^ ~ 1 done : 



n\ sin 2 i = ri2 cos r 



Et les solutions deviennent : 



E' 2ni cos i 

En m cos i + ?i2 cos r 

Eq n\ cos i — ri2 cos r 

E'0 ni cos i + ri2 cos r 

Le cas de la reflexion totale est obtenu pour cosr = 0, e'est-a-dire Eq = Eq. 
Pour calculer le transfert d'energie, nous pouvons utiliser la densite d'energie et 
le vecteur de Poynting pour un milieu : 



\(e-6 + b.h) 



U = 2 
S = Ex H 

La derivation de ces formules se fait de maniere similaire a celle pour le vide. 
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